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Abstract. Motivated by an earlier work on fractional-action cosmology with a periodic weight 
function 1 , we extend it by choosing a power-law weight function in the action. In this 
approach, we obtain a varying gravitational coupling constant. We then model dark energy 
in this paradigm and obtain relevant cosmological parameters. 



1. Dirac's Large Numbers Hypothesis and time running of G 

During our survey of various physical laws of nature, we came across a number of constants 
entangled with these laws. Historically it was Weyl [3] who initiated the idea of large numbers. 
However, it was Dirac who discovered an apparently unseen thread joining up these physical 
constants by a simple yet interesting law, viz., Law of Large Numbers. Using this law, Dirac 
arrived at his Large Number Hypothesis (LNH) which had profound influence on the world of 
physics. It generated a lot of activity both at theoretical and observational level with LNH as 
their starting point. Some works are centered around modification of Einstein's gravitational 



theory and related equations to adapt the G variation [U [5]. Some authors have also given 
counter-arguments against LNH for justifying and refuting that hypothesis characterize the 
second category of works in favor of it [6] and also testing the validity of it [7j. Even there is 
found a formal link between LNH and cosmological constant term A [8]. Weyl's large number is 
square root of Eddington's number N ~ 10 80 (which represents the total number of protons in 
the Universe) [9] and is supported by observational data [10J. Indeed Stewart showed that the 
ratio of the radius of the Universe and the electron is only two orders of magnitude smaller 10 
than that of Weyls number. There are some other large numbers like Jordan's number [TT] and 
the Shemi-zadeh's number [T2] which have also appeared in the literature. 

As originally stated by Dirac there are three dimensionless numbers in nature which can be 
constructed from the atomic and cosmological data: 

(i) The ratio of the electric to the gravitational force between an electron and a proton 
j^— ~ 7 x 10 39 . 

&m e m p 

(ii) The mass of that part of the Universe that is receding from us with a velocity v < c/2 
expressed in units of the proton mass is 10 78 . 
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(iii) The age of the Universe t, expressed in terms of a unit of time e 3 provided by atomic 
constants, is inversely proportional to Newton's gravitational constant. 

The above three statements imply in respective order 



e 2 



Gm e m p 



oc t, 



M oc i 2 , (2) 
Gott- 1 . (3) 

There is a serious problem between (3) and GR [H H3]: Einstein's theory requires G to be 
a constant. As was noted by Dirac, this inconsistency might be resolved if", we assume 
that Einstein theory is valid in a different system of units from those provided by the atomic 
constants. " Consequence of Dirac's LNH is the existence of a variable G cosmology. A beautiful 
account of the constants of nature appears in |14j . 



2. Framework 

Recently Nabulsi |15| examined the Friedmann-Robertson- Walker (FRW) spacetime in the 
context of fractional actiondike variational approach or fractionally differentiated Lagrangian 
function depending on a fractality parameter \ to model nonconservative and weak dissipative 
classical and quantum dynamical systems. 

We apply the technical machinery of differential geometry to general relativity in the context 
of fractional actiondike variational approach. Following Einstein's arguments, we start by 
representing the gravitational field by an affine connection on a curved manifold and assume 
that the free fall corresponds to geodesic motion. Thus the classical Lagrangian of a curve is 
defined by L = (l/2)gijx l x 1 , where x 1 = dx % jdt and is the metric tensor. This gives the 
action 

T 



where m is the particle's mass. Here n(x, t) is the generalized weight function and if we choose 
it to be a periodic function, we recover the results in pp. This weight function can produce a 



friction force in classical dynamics for the test particle and as we will show, it also yields a time 
varying G. The modified resultant fractional Euler-Lagrange equations are 

x l + dt log(/x(x, t))± l + T l ij ± i ± j = 0. (5) 

We follow the procedure of [15] and assume that the test particle's velocity is approximately 
zero. Furthermore, the metric and its derivative are approximately static. Applying these 
assumptions to the spatial components of (5) gives: 



dii* + dtfxdiv* + R tt = 0. (6) 

In the above equation v l = dx l /dt and 7* = dv l /dt are the particle's velocity and acceleration 
respectively and drf = —AitGp, which is Poisson's equation in Newtonian gravity. Here p is 
the mass density in gravitational interaction and G is Newton's gravitational constant. The 
resultant modified gravitational constant is 



G' = Gil - fti«g(M(x,t)) qA = G Aa (7) 

V AtrGp > 

If n(x,t) = t then (7) describes Dirac's model [T3]. Thus the change in G in a spatially flat 
(D + 2)-dimensional FRW Universe is 

AG = D(D + l)Hd t \og(p(x,t)) ^ 

8-7T/9 

where H = c^(log(a(i)) is the Hubble parameter. The Friedmann equation for such a model 
with metric gij = diag(l, —a 2 (t)rjAB), A, B = 2, D + 2, with t]ab = diag(l, 1) is 

D ( D 2 +1) H* = 8,G P + A. (8) 

One thinks that if the gravitational constant changes slightly, as we have expressed it in this 
work, the Friedmann's equation must be replaced by an effective G. A simple example of effective 
G is the f(R) gravity. An interesting feature of these theories is the fact that the gravitational 
constant is both time and scale dependent. After some calculations, one can obtain a Poisson 
equation in the Fourier space and attribute the extra terms as effective gravitational constant 
G e ff [IS]- But now, we are working using an Einstein-Hilbert action therefore there is no need 
to insert effective G on the right hand side of the Friedmann's equation (8). 

Using the Bianchi identity and assuming both energy density p and the cosmological constant 
A as functions of time related by A(i) = Aop(t) n . Consequently we have 

where A n = A ( 87r ) n/(n ~ 1} and p n = ( 8n G ^- n) )V("-i). Substituting (9) in (8) we 
obtain 

# = - = c „Mx,*r /(2(1 - n)) , (10) 



where c n = J D ^ D \i^ and b n = 8irGp n + A n and < n < 1. By solving (10), we find 
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Figure 1. Variation of scale factor (13) for different values of n. 



a(t) = a exp (c n J M ( X , i) n/(2(1 " n)) rft) • 
For the remaining part of the paper, we use power law weight function as 

^( x ,t)= X t p , p>l. 



(11) 



(12) 



One can treat (12) as a generalized Dirac's hypothesis. Note that p = 1 leads to an inverse 
power of time relation for the gravitational constant i.e. G = ^f-. Substituting (12) in (11) 
yields 



where 



a(t) = ao exp(a# 



a = x W^ ^-n) t ^^2 + (p-2)n 



(13) 



(14) 



2 + (p-2)n' " " 2(1 -n) 

It is useful to discuss the behavior of the scale factor diagrammatically. For simplicity we take 
arj = 1 and a = — 1, and focuss only on Dirac's model i.e. G(t) = j. Figure 1 shows the 
behavior of the scale factor as a function of time for different values of n. 
The Hubble parameter is 



H = apt?- 1 = a p( - - ln(l + z 



1-- 



(15) 



Here redshift z is defined by 1 + z = . Figure 2 shows the behavior of the Hubble parameter 
as a function of time for a = 1. 




Figure 2. Hubble parameter against t for selected values of n. 



We note that the scale factor for n = | , a = | is o oc exp(y). Consequently, the scale factor 
would increase very rapidly for the positive branch, while it would start from infinity for the 
negative branch. Also the choice n = | , a = gives a oc exp(^p). This indicates that the 
Universe is evolving with a Hubble parameter H ~ —t and a scale factor a oc exp( z ^), which 
is close to singularity related to the origin and the fate of the Universe. The same phenomena 
occurs in the context of a new model of dark energy, the Quintom model |17|. 



2.1. Deceleration parameter 
We know that 

H , , 

« = - 1 -H2> ( 16 ) 

is the deceleration parameter. For an expanding Universe, q > — 1. Using (14) we can write 
(16) as 

1-/3 1 . . 

q = ~ 1 ~~f~Ul + £)' (17) 



For Dirac's hypothesis (p = 1), we have 



72 [ 

" 1+ 2-nln(l + z)- (18) 



Thus we conclude that the Universe can remain in an accelerating expansion for some values of 
z. The next diagram (Figure 3) represents the graph of q versus z for different values of n. It 
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Figure 3. Deceleration parameter q against redshift z for some values of n. 



shows that when n increases the deceleration parameter q recedes from the crossing line q = — 1 
0. Thus we can model cosmic acceleration in the present framework. 

2.2. The statefinder parameters for dark energy 

The statefinder parameters {r, s} were introduced to characterize primarily the spatially flat 
Universe (k = 0) models with cold dark matter (dust) and dark energy [18]. They were defined 
as 



(/3 2 -3/3 + 2) i" 2 ^ + /3 ((3/3-3)^ + a(3)a 



iH 3 



a 2 (3 2 



(19) 



r _l 2 (Z 3 - 1 ) (t" /3 /?-2t- /3 + 3a/3) 

3(g- i) ~ ~3 a/3 (30^/3 + 2/3-2) 



(20) 



From (19) and (20) we obtain the relation between r and s as 

_ a 2 ' 13 (2 + 7s/3y-4/3-9s + 16s/3 + 2 s 2 /3 2 - 13 s 2 /3 - 7 s/3 2 - 7 - 2 /3 y + 2 /3 2 + 11 s 2 + 2 y) 



(-/3 + 1 - s/3 + s + V y 



where 



y = v //32_ 2/ 3 + 2s/3 + l-2s + s 2 /3 2 -2s 2 /3 + s 2 . 

1 remembering that the < n < 1 



(22) 




Figure 4. The statefinder parameters for some values of n. 



We know that the cosmological constant possesses a fixed equation of state parameter 
(w = —1) and a fixed Newton's gravitational constant, hence {1,0} corresponds to ACDM. 
Moreover {1,1} represents the standard cold dark matter model containing no radiation while 
Einstein's static Universe corresponds to {—00,00} [T9] , 

Note that for Dirac's hypothesis p = 1, n = 0, = 1, we have {r, s} = {1, 0} which represents 
a static cosmological constant with cold dark matter. In Figure-4, we draw the pair {r, s} for 
different values of (5 = 2 (i-n) > a = ^ ^ s we can see ' ^ ne fi xe d gravitational constant occurs 
several times for some values of n. Also our model describes the standard cold dark matter and 
a static Einstein Universe for some values of n. 



2.3. Lookback time 

If a photon is emitted by a source at the instant t and received at time to, then the photon 
travel time or the lookback time t — to is defined by 

f z dz 

t-to = - , (23) 

Jo H(z)(l + z) 

where ao is the present value of the scale factor of the Universe. It is appropriate to calculate 
this time as a function of z. After substituting Eq. (15) in Eq. (23) we have 

t-t = - K 1 



a/3 Jo (log(l + z)) 1 - 1 //3(l + z)' 
Now we change the variable as x = log(l + z). Calculating the integral above, we find 

t-t = (- l ^±^)k (24) 
a 



The behavior of the lookback time as a function of the redshift for 



values q=-l,0<n<l 
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Figure 5. The behavior of the lookback time as a function of redshift for a = — 1, < n < 1. 



Figure-5 shows the behavior of the lookback time as a function of redshift for values a = — 1, < 
n < 1. 

2.4. Proper distance 

If a photon is emitted by a source and is received by an observer at time to then the proper 
distance between them is defined by 

f a " da f z dz 

d = ao L ^p = Jo irwi+zy (25) 

Using Eq. (15) in (25), we get 

Figure 6 shows the behavior of proper distance in the present model. 

2.5. Luminosity distance 

If L is the total energy emitted by the source per unit time and £ is the apparent luminosity of 
the object, then the luminosity distance di = (4^) 2 evolves as 

,/-, n z + l ( ln(l + z)\-i+l 

Figure (7) shows the variation of the luminosity distance as a function of the exponent 
< n < 1. 
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Figure 6. The behavior of the proper distance as a function of redshift for a = — 1, < n < 1. 




Figure 7. The behavior of the luminosity distance as a function of redshift for a = —1,0 < 
n < 1. 




Figure 8. The behavior of the angular diameter as a function of redshift for a = — 1, < n < 1. 



2.6. Angular diameter 

The angular diameter distance cIa is 

Figure-8 shows the variation of the angular diameter as a function of the exponent < n < 1 . 



2. 7. State parameter of dark energy 

In this subsection we discuss the dynamics of the dark energy. We know that in a flat FRW 
model with a perfect fluid of pressure p and energy density p, the barotropic index w(t) is defined 
by 

-(«)— <W 
which can be can be derived from the continuity equation. Using (9), (10) with (29) we obtain: 

' n pn+2— 2 n 

W (t) = -1 ; r npx 2{ - 1+n) t 2(-l+n) . (30) 

3(-l + n)c„ 

For simplicity we take D = 2, x << 1 an d G = Like the previous sections, we are 
interested in expressing w as a function of z. We substitute /3 = 2 ^~^" and from Eqs. (24) 
and (30) we have 

1 ^^/r-loga + z)!^^^ 



3 (-1 + n) c n 



a 



The behavior of the barotropic index w vs. z 



w(z) 




Figure 9. The behavior of the barotopic index w vs. z 



Since ^ x P 2(i-4" = ^ therefore the previous equation reads 



w(t) = -1 + 



3a (—1 + n) c 

Now using the values of p n , b n , c n and A n , we can write 

w(z) = -1 - d n \n(l + z) , 

where 



n PX 2(1_n) log(l + z), 



dn = -7=( 



2p n n x 



n,2—n 



(31) 



(32) 



v/3 v Ao(l-n)' 

Also we observe from the previous figures, the model parameters do not depend on p 
significantly and since the fractal parameter x << 1 therefore we limited ourselves only to 
p = 1. From (30) it is obvious that the barotropic index crosses the dark energy limit w(t) = —1 
for current era z = several times for different values of x an d n. For instance, using typical 
values x = 0.05, n = 0.5, A = .01826, (30) becomes 

w {t) = -1 - In (1 + z) . (33) 

Figure-9 shows the behavior of the barotopic index w against z. The pressure may be written 



as 



where 



P(z) = -Pn(ln(l + z) + d„ln(l + z) — ), 

Pn=Pn{ ) 

V a' 



(34) 
(35) 



2.8. Consistency of the model G{t) = % with observational constraints 

From observational data, as it was shown in [20], there is a secular decrease of the gravitational 
constant G(t). This time rate of change is of order ^ = (—5.9 ± 4.4) x 10~ 14 yr~ 1 which 
was obtained from planetary data analysis. Thus, in our treatment on fractional cosmology 
focusing on the model G(t) = j, we can accept it as a reasonable model in good agreement with 
observational constraints |21j . 

Final remarks 

This work was motivated by a more recent work on fractional-action cosmology with a periodic 
weight function in the action pQ which we extended using a power-law weight function. By doing 
so we obtained a varying gravitational coupling constant and also a time varying cosmological 
constant. For both of these 'constants', we used different ansatz to analyze the behavior 
of other parameters. We also modeled dark energy in this paradigm and obtained relevant 
cosmological parameters including distance parameters and statefinder parameters. The former 
helps in determining the nature of dark energy while the latter one classify different dark energy 
candidates. 
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